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, The constrained Hamiltonian systems admitting no gauge conditions are considered. 

The methods to deal with such systems are discussed and developed. As a concrete 
application, the relationship between the Dirac and reduced phase space quantizations 
is investigated for spin models belonging to the class of systems under consideration. 
It is traced out that the two quantization methods may give similar, or essentially 
' different physical results, and, moreover, a class of constrained systems, which can be 

quantized only by the Dirac method, is discussed. A possible interpretation of the 
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1 Introduction 



Nowadays there are two main methods to quantize the Hamiltonian systems with first class 
constraints: the Dirac quantization [0] and the reduced phase space quantization 0, while 
the path integral method 0, |[ and the BRST quantization, being the most popular method 
for the covariant quantization of gauge-invariant systems |Q, are based on and proceed from 
them 0, ^, ^. The basic idea of the Dirac method consists in imposing quantum mechanically 
the first class constraints as operator conditions on the states for singling out the physical 
ones The reduced phase space quantization first identifies the physical degrees of freedom 
at the classical level by the factorization of the constraint surface with respect to the action 
of the gauge group, generated by the constraints. Then the resulting Hamiltonian system 
is quantized as a usual unconstrained system 0|. Naturally, the problem of the relationship 
of these two methods arises. It was discussed in different contexts in literature (see [[7[-p3|, 
and references therein), and it was demonstrated for the concrete examples that these two 
methods give, in general, different results. The corresponding quantum systems may differ, 
e.g., in their energy level structure, or in the eigenvectors of their Hamiltonians 0, |10 



Nevertheless, there is an opinion that the differences between the two quantization methods 
can be traced out to a choice of factor ordering in the construction of various physical 
operators. 

In the present paper we investigate the relationship of the two methods of quantization 
for the special class of Hamiltonian systems with first class constraints corresponding to dif- 
ferent physical models of spinning particles. The specific general property of the considered 
examples of constrained systems is the following: their constraints generate SO (2) transfor- 
mations and, hence, corresponding gauge orbits topologically are one-spheres S^. This fact 
implies that these systems do not admit gauge conditions. Nevertheless, for the first three 
systems we can construct the corresponding reduced phase spaces. We shall demonstrate 
that the two methods of quantization give coinciding physical results only for the first of 
these systems, but they give, in general, essentially different physical results for two other 
systems. 

Since in the present paper we consider the constrained systems which do not admit gauge 
conditions, we should use a general geometrical approach to the Dirac-Bergmann theory of 
the constrained systems for the construction of the reduced phase spaces. This approach is 
reviewed in Section 2 of our paper. Note that from the point of view of this approach, it 
is not necessary to use the notion of the weak equality of the functions on the phase space, 
and this notion is not used in our paper at all. 

The first concrete example, we consider, is the plane spin model, which is a subsystem of 
the (3-1-1) -dimensional models of massless particles with arbitrary helicity fl^, and of the 
(2-1-1) -dimensional relativistic models of fractional spin particles |]15|. This model is related 



also to the simplest nonrelativistic model of the anyon [|16]. The classical and quantum 
treatment of the model having a cylinder S*^ x R as the initial phase space, is performed 
in Section 3. Here we trace out an interesting analogy in interpretation of the situation 
with nonexistence of a global gauge condition for this simple constrained system with the 
situation taking place for the non-Abelian gauge theories fl^ . 



The second model, which is considered in Section 4, is the rotator spin model, being a 
subsystem of the (3-1-1) -dimensional model of the relativistic massive particle with integer 
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spin [1^. The nontriviality of the model reveals itself in the topology of the constraint 
surface diffeomorphic to the group manifold of the Lie group SO (3). 

Section 5 is devoted to the consideration of the top spin model, which, in turn, is a 
subsystem in the non-Grassmannian model of the relativistic massive particle with arbitrary 



integer or half-integer spin |T9|. This system is also topologically nontrivial and has the 
cotangent bundle T*S0(3) as the initial phase space. 

In section 6 we summarize the results, present conclusions on the relationship of the two 
methods of quantization and formulate a possible interpretation of the sense of gauge degrees 
of freedom for constrained systems. Here we also discuss a broad class of pseudoclassical 



models p0[-pj], i.e., models containing Grassmann variables, whose general feature consists 
in the presence of the constraints, nonlinear in Grassmann variables and having no nontrivial 
projections on the unit of Grassmann algebra. For such systems it is not possible to construct 
the reduced phase space in principle, and they can be quantized only by the Dime method. 

Appendixes A and B are devoted to the presentation of the method of dependent coor- 
dinates used in the paper. Actually this method is, in a sense, a reformulation of the Dirac 
bracket method, but from our point of view it has, in comparison with the former, greater 
mathematical definiteness, and turns out to be more convenient for our purposes. 

Appendix C collects the facts about three dimensional rotation group, which are needed 
in the main text. 

Everywhere in the text repeated indices imply the corresponding summation. 



2 Geometrical Background of the Dirac— Bergmann 
Theory 

As it has been noted in the Introduction, the systems we consider in the present paper do 
not admit gauge conditions. To construct the reduced phase space for them we should use a 



general geometrical approach (see [^, and references therein). This approach together with 
its relation to the usual Dirac-Bergmann formalism is briefly described below. 

Let (M , u'^) be an m-dimensional symplectic manifold, and be an n-dimensional 
submanifold of M. Let l : N ^ M he the inclusion mapping. The two-form a;^ = i*u^ is 
obviously closed. Suppose that has a constant rank, then the set 

N 



S^N = {xeTN \ i{x)uj^ = 0} 



is a distribution on A^, called the characteristic distribution of cu . As cu is closed, £^n 
is integrable and by the Frobenius theorem defines a foliation jF^jv on A^. Identifying all 
points of a leaf we get the quotient space P = N/J^^n. Assume that P is a manifold with 
the canonical projection tc : N —>■ P being a submersion. In other words, we suppose that 
TT : A^ ^ P is a fibered manifold. In this case it can be shown that there exists a closed 
nondegenerate two-form u^, such that 

7r*o;^ = cu^. (2.1) 
Thus, the manifold P has a natural structure of a symplectic manifold. 
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In the Dirac-Bergmann theory the submanifold is a manifold defined by equations 

ipa = 0, a = 1, . . . ,m — n. 

It means that 

N = {p e M \ tpaip) = 0, a = 1, . . . ,m - n}, 

where ijja are differentiable functions on M, such that the mapping ip : M ^ R™~", defined 
by p — ^ (V'i(p); • • • ; '^m-n{p)) IS of rank m — n for every p & N. 
Let X G TpN, then 

d4'a{i^*{x)) = 0, a = 1, . . . ,m — n. 
On the other hand, suppose that y G T^(^p)M, p E N, and 

#a(z/) = 0, a = l,...,m-n, (2.2) 

then there exists a unique vector x G T^A^, such that y = l^{x). 

Recall that the Hamiltonian vector field Xf on M, corresponding to the function / on 
M, is defined by 

t{Xf)u;^' = df, (2.3) 
and for the Poisson bracket of two functions, / and g, on M we have the expression 

{/, g} = -cu'^Xf, X,) = Xf{g) = -X,(/). (2.4) 

Introduce the matrix valued function A = ||Aaf,||, where 

and define 

Suppose, that the matrix A is degenerate at some point p E N . Hence we can find a 
nontrivial set of real numbers c°, a = 1, . . . ,m — n, such that 

c'^Aabip) = 0. 

From Eq. (p. 41) we get 



c''Xa,(p){A) = dlpu(p){.c''Xa) = 0, 

where we have denoted X^^ through Xa- Thus, there exists a vector x G TpN, such that 

i,{x) = c'^Xa^^p). (2.5) 

Let us show that the vector x belongs to S^n. Indeed, for any x' G TnN we have 

[t{x)uJ^]{x')=Uj''{x,x')=UJ^'mx),L,{x')). 

Using Eq. (|2.5|), we get 

[t{x)uj'']{x') = C^CO^'iXayi^x')) = ediJa{i.{x)) = 0. 
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That was to be demonstrated. 

On the other hand, suppose that x e Ei^n. li x E TpN, then for any x' e TpN we have 

0=[z{x)u'']{x') =u''mx),L,{x')) = [^(t,(a;)V^](t,(a;')). 

Note that if a one-form /j, e T*{p)M is such that 

Ii{l^{x)) = 

for all X e TpN, then 

for some constants c", a = 1, . . . ,m — n. Hence, we have 
and 

Thus, we see that the elements of E^n and the null vectors of the matrices A(p), p E N, are 
in one-to-one correspondence, and the two-form cu^ is of constant rank if and only if the 
rank of A(p) is independent of p. 

Recall that in the Dirac-Bergmann theory the functions ipa are constraints, the subman- 
ifold N is called the constraint surface, and the symplectic manifold P is called the reduced 
phase space. 

A function / on M is said to be of first class, if 

{/, ^a} = fa^b 

for some functions A constraint being a function of first class, is called the first class 
constraint. Suppose that the matrix A{p) is of rank k at any point p G M. Let the indices a, 
P, 7 take values from 1 to k, while the index /j, takes values from k+1 to m—n. Without a loss 
of generality, we can consider, in fact locally, that the submatrix ||Aq,^(p)|| is nondegenerate 
at any point p E M. In this case we can introduce the functions A"^ satisfying the relation 

A"\p)AM=d-^ 

for all p e A^. It is clear that the functions 

</';. = V'm-{^M>V'JA"% 

form an equivalent set of constraints defining the constraint surface A^. 

The constraints 0^ are first class constraints. It can be easily shown that the Hamiltonian 
vector fields X^^ are tangent to A^, and these fields restricted to A^ form a basis of the 
characteristic distribution £^n. The main property of the constraints which certainly 
follows from their definition, is that the matrix formed by the Poisson brackets of 0q, is 
nondegenerate at any point of the constraint surface. The constraints (pa are called second 
class constraints. 
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If the matrix A(p) is of zero rank at any point p & N, we have only first class constraints. 
In this case 

{^a,i^b}=^crab (2.6) 

for some functions f^ab on M. If the matrix A{p) is nondegenerate at any point of N we deal 
with second class constraints. In this case the corresponding characteristic distribution is 
trivial, and the reduced phase space coincides with the manifold A^. A convenient method to 
work with second class constraints is the Dirac bracket method, or the method of dependent 
coordinates (see Appendix B). Note that in this case it is convenient to use upper indices 
for the constraints. 

Let us discuss now the procedure of constructing the reduced phase space for the case of 
first class constraints, according to the Dirac-Bergmann theory. 

Below we identify the tangent space TpN with the set of vectors from T^{jp)M, satisfying 
Eq. (p.2|) , and consider the submanifold A^ simply as a subset of M. If the matrix A(p) 
is of rank zero at any point p & N, then the Hamiltonian vector fields Xa belong to the 
characteristic distribution S^^n, and generate it. In the case when the functions f^ab from 
Eq. ( |2.(j| ) are constant on M, we can try to find a group of canonical transformations G such 
that the functions ipa be its generators. This group is called the gauge group, generated 
by the constraints ipa- If such a group exists, the leaves of the foliation, defined by the 
characteristic distribution, are its orbits. 

Suppose that the fibered manifold tt : N P has a global section s : P —>■ N. The 
differentiable mapping s by definition satisfies the relation 

TT o s = idp. (2.7) 

It is clear that the image s{P) of the section s can be considered as a submanifold of A^ that 
is diffeomorphic to P. Furthermore, from Eqs. (|2.1| ) and ( p.7| ) it follows that 

* N P 
S OJ = OJ . 

Hence, 

where uj^^^'^ is the restriction of the two-form oj^ to the submanifold s(P), and s is considered 
as a mapping from P to s{P). Thus the submanifold s{P) is a symplectic manifold that is 
symplectomorphic to P, and we can identify P with s{P). Choosing different global sections 
of the fibered manifold tc : N —>■ P, we come to different realizations of the reduced phase 
space as a submanifold of A^. Note that A^ is a submanifold of the symplectic manifold M, 
hence we actually have realizations of the reduced phase space as a submanifold of M. 

Assume now that the reduced phase space, realized as a submanifold of M , is defined by 
the equations 

i'a = 0, x" = 0, a = l,...,m-n. 

The functions x° are called gauges, or gauge conditions. They can be considered as additional 
constraints. Denote the full set of the constraints by ^z^", a = 1, . . . , 2(m — n): 
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As P is a symplectic submanifold of M, the matrix with the matrix elements = 
L*{ip"',ilj'^}{p) is nondegenerate at any point p & P. From this it follows that the matrix 

Aa'(p) = ^*{^a,/}(P) (2.8) 

must be of rank m — n at any point oi p E P. 

Note that a section s allowing to identify the reduced phase space P with a symplectic 
submanifold s{P) of the initial symplectic manifold M may not exist. In this case we say 
that the system under consideration does not admit gauge conditions. In such a situation we 
can always consider a set of local sections Si : Ui ^ N, covering P. For any i the set Si{Ui) 
is symplectomorphic to the open symplectic submanifold Ut. The total reduced phase space 
can be considered as the result of natural patching of submanifolds Sj(f/j). Actually it is 
more convenient to construct the reduced phase space by the direct geometrical construction 
described above. It is the method we use in the present paper. 

Note also that even if a global section of the fibered manifold it : N ^ P exists we 
cannot, in general, consider s{P) as a submanifold of M, defined by equations. We shall 
encounter such a situation in the next section. 



3 Plane Spin Model 



The first model we are going to consider is the plane spin model 11^ , [15| . This is a very 
simple model, but nevertheless here we come across peculiarities relevant for a more general 
case. Therefore, we give for this model a more detailed consideration than it deserves from 
a first sight. The initial phase space of the model is a cotangent bundle T*S^ of the one- 
dimensional sphere S*^, that is a cylinder S"^ x R. It can be described locally by an angular 
variable < ip < 2tt and the conjugate momentum S" G R. Here the angle variable cannot be 
considered as a global coordinate. The symplectic two-form u in terms of the local variables 
<f, S has the form 

u = dS A dip. 

Thus, we have locally 

{V,S} = 1. 

Actually, any 27r-periodical function of the variable (p that is considered as a variable, taking 
values in R, can be considered as a function on the phase space, i.e., as an observable, and 
any observable is connected with the corresponding 27r-periodical function. Therefore, we 
can introduce the functions 

gi = cosy?, q2 = sirup (3.1) 
as the functions on the phase space of the system. For these functions we have 



{qi,q2} = 0, 
{qi,S} = -q2, {q2,S} = qi. 

The functions qi and q2 are not independent, and satisfy the relation 

qf + ql = 1. 



(3.2) 
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Note, that any function on the phase space can be considered as a function of dependent 
coordinates qi, q2 and 5*. These coordinates will be taken below as the quantities, forming a 
restricted set of observables whose quantum analogs have the commutators which are in the 
direct correspondence with their Poisson brackets. 

We come to the plane spin model by introducing the 'spin' constraint 

-ilj = S-e = 0, (3.3) 

where 6 is an arbitrary real constant. 

Consider now the Dirac quantization of the system. To this end let us take as the Hilbert 
space the space of complex 27r-periodical functions of the variable (p with the scalar product 

(<l>i,$2) = — / $i((/p)$2(¥^)rf</^. (3.4) 
Ztc Jo 

The operators qi and q2, corresponding to the functions qi and q2, are the operators of 
multiplication by the functions cosip and siny?, respectively, i.e.. 

The operator S is defined by 

where c is an arbitrary real constant. It is clear that the operators qi, q2 and S are Hermitian 
operators with respect to the scalar product, defined by Eq. ( p.4|) . 
One can easily show that the relations 

[^1,^2] = 0, 
[qi,S] = -iq2, [q2,S]=iqi 

are valid for any value of c in agreement with Eq. ( p.2|) . 

The quantum analog of the constraint ( |3.3| ) gives the equation for the physical state wave 
functions: 

iS-9)%hys = 0. (3.5) 
Decomposing the function ^phys{<p) over the orthonormal basis, formed by the functions e*'^'^: 

fcez 

we find that equation (|3.5| ) has a nontrivial solution only when 

c = 9 + n, 

where n is some fixed integer, n G Z. In this case the corresponding physical normalized 
wave function is 

%hysip) = e'"'^. (3.6) 

Here the only physical operator 0], i.e. an operator commuting with the quantum constraint 
ip, is S, and it is reduced to the constant 9 on the physical subspace 
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Now we come back to the classical theory in order to construct the reduced phase space 
of the model. Let us show that for the surface, defined by Eq. (|3.3| ), there is no 'good' gauge 
condition, but, nevertheless, the reduced phase space of the system can be constructed. In- 
deed, it is clear that the one-parameter group of transformations, generated by the constraint 
ip, consists of the rotations of the phase space. This group acts transitively on the constraint 
surface, and we have only one gauge orbit, which is the constraint surface itself. The gauge 
conditions must single out one point of an orbit. In our case we have to define only one 
gauge condition, let us denote it by x- The function x must be such that the equations 

^ = 0, X = 

determine a set, consisting of only one point, and in this point we should have 

(see Section 2). Recall that any function on the phase space of the system can be considered 
as a function of the variables ip and S, which is 27r-periodical with respect to Thus, we 
require the 27r-periodical function x(v^? 'S')|5=e turns into zero at only one point = ipQ from 
the interval < ^9 < 27r. Moreover, we should have 



7^0. 



It is clear that such a function does not exist. Nevertheless, we have here the reduced phase 
space that consists of only one point. Therefore, the reduced space quantization is trivial: 
physical operator S takes here constant value 9 in correspondence with the results obtained 
by the Dirac quantization method. When the described planar spin model is a subsystem of 
some other system, the reduction means simply that the cylinder T*S^ is factorized into a 
point, where S = 6, and that wave functions do not depend on the variable if. 

Let us point out one interesting analogy in interpretation of the situation with nonexis- 
tence of a global gauge condition. Here the condition of 27r-periodicity can be considered as 
a 'boundary' condition. If for a moment we forget about it, we can take as a gauge function 
any monotonic function xi'^^ S), x ^ such that xlv^o? 6^) = at some point ip = ipo, and, 
in particular, we can choose the function x(v^) S) = ip. The 'boundary' condition excludes 
all such global gauge conditions. In this sense the situation is similar to the situation in the 
non-Abelian gauge theories where without taking into account the boundary conditions for 
the fields it is also possible to find global gauge conditions, whereas the account of those 
leads, in the end, to the nonexistence of global gauge conditions (see, e.g., [0). 

Now we are going to describe the system using dependent coordinates (see Appendix B). 
To this end let us take two-dimensional vectors q and p, where p is defined as 

Pi = —Ssimp, p2 = Scosip. 

Let eij be a skew symmetric tensor, normalized by the condition = 1. For any two 
dimensional vectors a and b we use the notation 

a X b = eijaibj. 
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It is clear that 

-S" = eijQiPj = qxp. 
The variables q and p satisfy the relations 

q2 = 1^ qp = 0, 

and can be considered as dependent coordinates in the phase space. For nonzero Poisson 
brackets we have the following expressions: 

{Pi,Pj} = PiQj - QiPj, {(li,Pj} = Sij - mj- (3.7) 



Hence, the Hamiltonian vector field, corresponding to constraint ( |3.3| ) has, according to 
Eq. (ra), the form 

^ d d 

oq op 

This vector field generates the characteristic distribution of the reduction of the symplectic 
two-form to the constraint surface, given by Eq. ( p.3|) . The corresponding leaves are the 
orbits of the one-parameter group of canonical transformations, generated by the function 
ip. To find them we can use Eq. ( p.lOp , that leads to 

q(r) = q cos t + S^^p sin r, p(r) = p cos r — Sq sin r. 

It is clear that in accordance with the above consideration the constraint surface consists of 
only one leaf, and the reduced phase space is a point. 

Concluding this section, let us write the corresponding Lagrangian for the described 
Hamiltonian system: 

L = -eqxq-X{q^ -I), (3.8) 

where A is a Lagrange multiplier. It is not difficult to show that this Lagrangian leads to 
the Hamiltonian system with two second class constraints 

^i = qp = 0, 1P2 = \{q^ - 1) = 0, 
and one first class constraint 

%l) = qxp — 6 = Q. 

Using the notion of the Dirac bracket (see Appendix B), it is easy to get convinced that the 
Poisson brackets for the variables qi and p,, restricted to the submanifold, defined by the 
constraints ipi and %l)2i coincide with those, given by Eq. ( |3.7| ). It is convenient further to use 
the same letters for coordinates and their restrictions to the submanifold under consideration. 
The exact meaning of the coordinates is always clear from the context. 

Using parameterization ( |3.1| ), we come from Lagrangian ( |3.8[ ) to the Lagrangian 

L = Oip. (3.9) 

Note that Lagrangian (|3.9|) is present as a part of a complete Lagrangian in the nonrelativistic 
model of the anyon . 

Concluding this section, stress once more that for the considered simplest system the 
both quantization methods give coinciding results. 
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4 Rotator Spin Model 



In this section we consider the rotator spin model [|1^ . The initial phase space of the system 
is described by a spin three- vector S and a unit vector q, 



being orthogonal one to the other, 

qS = 0. 

The variables qi and St, i = 1,2,3, can be considered as dependent coordinates in the 
phase space of the system. The Poisson brackets for these coordinates are 

{liy Qj} = 0, {Si, Sj} = eijkSk, {Si, Qj} = eijkQk- 
Introducing the vector 

p = S X q, 

we get for its components the following Poisson brackets: 

{Pi,Pj} = PiQj - QiPj, {(li,Pj} = Sij - qiQj, 
that is of the same form as Eq. ( pj.T] ), and 

{Si,pj} = tijkPk- 
Note that the vector p satisfies the relation 

qp = 0. 

Actually, the variables and pi can be considered as another set of dependent coordinates 
in the phase space. 

Using the expressions for the Poisson brackets of the dependent variables, we find the 
following expression for the symplectic two-form: 



u = dpiA dqi = d{eijkSjqk) A dqi. 



(4.1) 



Note that the phase space under consideration is the internal phase space of the one-mode 
relativistic string model . 

Introducing the spherical angles ^p, ^ {0 < ip < 2tt, < < tt) and the corresponding 
momenta Pip,p-d G R, we can write the following parameterization for the vectors q and p: 



cos if sin 
sin sin d 

cos'd 



P 



cos V2 COS iip^Q : — -p^ ^ 

smw 

n COS 09 

sm vj cos Vp^ + p^p 



sint?' 



V 



— sin dp^ 



(4.2) 



and for the symplectic two-form we get the expression 

ijj = dp^ A d-d + dp^ A dip. 
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From this relation we conclude that the initial phase space of the system is symplectomorphic 
to the cotangent bundle T*S'^ of the two-dimensional sphere S"^, furnished with the canonical 
symplectic structure. 

The rotator spin model is obtained from the initial phase space by imposing the constraint 

^ = ^{S'-p') = 0, p>0, (4.3) 
fixing the spin of the system. This can be obtained also starting from the Lagrangian 

L = _p^-A(g2-1), 

where A is a Lagrange multiplier. Indeed, the Hamiltonian description of the system with 
this Lagrangian results in two second class constraints 

^1 = = 0, ^2 = liq' - 1) = 0, 



and a first class constraint (|4.3| ) with S = q x p. The reduction to the surface, defined by 
the second class constraints ipi and ip2, leads to the Poisson brackets which coincide with 
those given above. 

Using the Dirac method, we quantize the model in the following way. The state space is a 
space of the square integrable functions on the two-dimensional sphere. The scalar product 
is 

($i,$2)=/ <l>i(v9,^9)<l>2(<^,?9) smMM(p. 
Parameterization ([4.2|) allows us to use as the operator S the usual orbital angular momen- 



tum operator expressed via spherical angles The wave functions as the functions on a 
sphere are decomposable over the complete set of the spherical harmonics: 

oo I 
1=0 m=-l 

and, therefore, the quantum analog of the first class constraint ( [4.3| ), 

(52-p2)$^^^^ = 0, (4.4) 
leads to the quantization condition for the constant p: 

p'^ = n{n + l), (4.5) 



where > is an integer. Only in this case equation ( f4.4| ) has a nontrivial solution of the 
form 

m=—n 

i.e., with the choice of ( [4.5|) we get the states with spin equal to n: 
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Thus, we conclude that the Dirac quantization leads to the quantization (O) of the param- 
eter p and, as a result, the quantum system describes the states with integer spin n. 

Let us turn now to the construction of the reduced phase space of the system. Then we 
perform the reduced phase space quantization and reveal that this method of quantization 
gives results physically different from those obtained by the Dirac quantization method. 

The constraint surface of the model can be considered as a set composed of the points 
specified by two orthogonal normed three-vectors. Each pair of such vectors can be sup- 
plemented by a unique third three-vector, defined in such a way that we get an oriented 
orthonormal basis in three dimensional vector space. It is well known that the set of all 
oriented orthonormal bases in three dimensional space can be smoothly parameterized by 
the elements of the Lie group SO (3) (see Appendix C). Thus, the constraint surface in our 
case is diffeomorphic to the group manifold of the Lie group SO (3). 

As in the preceding Section, to find the one parameter group of canonical transformations, 
generated by the constraint if), we can use Eq. ( [B.10| ). That gives the result 

q{T) = q cos(Sr) + ^Sxq sin(Sr), 5(r) = 5, S= V^. 

Hence, we see that the transformations we are looking for, are the rotations about the di- 
rection, given by the spin vector. Thus, in the case of a general position the orbits of the 
one-parameter group of transformations under consideration are one dimensional spheres. 
Note, that only the orbits, belonging to the constraint surface where S = p 0, are in- 
teresting to us. It is clear that an orbit is uniquely specified by the direction of the spin 
three-vector S whose length is fixed by the constraint tp. As a result of our considera- 
tion, we conclude that the reduced phase space of the rotator spin model is the coset space 



SO(3)/SO(2), which is diffeomorphic to the two-dimensional sphere S"^ |28]. Due to the 
reasons discussed in the preceding section there is no gauge condition in this case either. In 
fact, since S0(3) is a nontrivial fiber bundle over S"^, we can neither find a mapping from S"^ 
to SO (3) whose image would be diffeomorphic to the reduced phase space. In other words, 
in this case the reduced phase space cannot be considered as a submanifold of the constraint 
surface. 

Our next goal is to write an expression for the symplectic two-form on the reduced 
phase space. We shall give here two considerations of this question, using independent and 
dependent coordinates. 

Convenient independent coordinates in the initial phase space can be introduced here 
with the help of the relations 

cos 7 cos COS (p — sin 7 sin cp 
cos 7 cos sin cp + sin 7 cos p 
— cos 7 sin 

S < 00. In terms of these coordinates the 
symplectic two-form u takes the form 

uj = d{Sd-f + S cos Mp). (4.6) 

From this relation we get the following expressions for the nonzero Poisson brackets 

{l,S} = l, 

cos??} = S*^"*^, {cost?, 7} = S^"*^ cos-i?. 





/ 


sin -(9 cos 99 \ 




s = s 




sin ^9 sin 99 , 






v 


cost? / 




where < < vr. 





< V?,7 < 27r, 


and < 
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The restriction of the symplectic two-form uj to the constraint surface, which we denote by 
the same letter, has the form 

UJ = pdicos'ddip). (4.7) 

It is clear that the transformations of the group, generated by tp, act on the independent 
coordinates as follows: 

S{t) = S, ■d{T)=d, (p{t) = (f, 7(r) = 7 + S'r. 

Thus, we see that the variables -d, ^ can be considered as coordinates in the reduced phase 
space, and the symplectic two-form on the reduced phase space is given by ([4.71). 

The same results, but in a more transparent form, are obtained using dependent coor- 
dinates. It is clear that we can consider the variables Si as dependent coordinates in the 
reduced phase space, thus the symplectic two-form on it may be expressed in terms of them. 
Note, to this end, that the vectors g, s = S/S and q x s form an orthonormal basis, and 
the differentials of the vectors q and s satisfy the relations 

qdq = 0, sds = 0, (4.8) 
sdq + qds = 0. (4.9) 

Write the restriction of the symplectic two-form 1|) to the constraint surface as 

u! = p{ds X q) A dq + p(s x dq) A dq. 

From Eq. ( [4.8| ) it follows that we can write 

dq = as + b{qx s), (4.10) 
ds = fq + g{qx s), (4.11) 

where a, b, f and g are some one-forms. From Eq. ( [4.10| ) we get 

s X dq = bq, 

and, therefore, 

(s X dq) Adq = 0. 

Further, from Eq. (^4.11|) we have 

dsx q = gs. (4.12) 

This equality implies 

{ds X q) Adq = f A g, 
where we have used Eq. ( [4.9| ). From Eq. ( ^.11| ) we have also 

(s X ds) Ads = -2/ A g, 

that gives finally 

uj = — ^(S X dS) AdS. (4.13) 
2p^ 
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Thus, we see that the dependent coordinates S\ 

= p\ (4.14) 

in the reduced phase space of the system provide a reahzation of the basis of the Lie algebra 
so(3): 

{Si,Sj} = eijkSk. (4.15) 

The quantization on the reduced phase space can be performed with the help of the 
geometric quantization method proceeding from the classical relations ( [4.13| )- (|4.15| ). This 
was done in detail earlier (see, e.g., Ref. [^), and we write here only the final results of this 
procedure. The constant p is quantized: 

p = j, 0<2jGZ, (4.16) 

i.e., it can take only integer or half- integer value, and the Hermitian operators, corresponding 
to the components of the spin vector, are realized in the form: 

where 

or, in terms of the dependent coordinates, 

_ Si - 182 
P + S3 ' 

Operators ( [4.17| ) act in the space of holomorphic functions f{z) with the scalar product 



in which the functions 



2j + l r r Mz)Mz) 



ij'^(xz^+"', m = -j,-j + l,...,j. 



form the set of eigenfunctions of the operator 5*3 with the eigenvalues S3 = m. These 
operators satisfy the relation 



instead of classical relation ( [4.14|) , and, therefore, we have the (2j+l)-dimensional irreducible 
representation Dj of the Lie group SU(2). 

Thus, we see that for the rotator spin model the reduced phase space quantization method 
leads to the states with integer or half-integer spin, depending on the choice of the quantized 
parameter p, and gives in general the results physically different from the results obtained 
with the help of the Dirac quantization method. Let us stress once again here that within 
the Dirac quantization method in this model the spin operator S has a nature of the orbital 
angular momentum operator, and it is this nature that does not allow spin to take half- 
integer values [^. In other words, in the Dirac quantization we quantize the system with 
the phase space being the cotangent bundle to the two-dimensional sphere, defined by the 
relation = 1, while in the reduced phase space quantization we start with the phase space 
being a two-dimensional sphere of the radius coinciding with the length of the vector S. It is 
this difference between the phase spaces that leads to the different results of quantizations. 
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5 Top Spin Model 



In this section we consider the top spin model [|T9l. The initial phase space of the model is 



described by the spin three-vector S, and by three vectors ej, forming a right orthonormal 
basis in R^: 

(see Appendix C). Denote the components of the vectors Cj with respect to the canonical 
basis of by Eij. The components Si of the vector S and the quantities Eij form a set of 
dependent coordinates in the phase space of the system. The corresponding Poisson brackets 
are 

{E,,,Em} = 0, (5.2) 
{Si, Ejk} = eikiEji, {Si,Sj} = eijkSk- (5.3) 



The vectors Cj, satisfying Eq. ( ^.1|) form a right orthonormal basis in R^. The set of all 
such bases can be identified with the three-dimensional rotation group (see Appendix C). 
Taking into account Eqs. ( p.2| ) and ( |5.3| ) we conclude that the initial phase space is actually 
the cotangent bundle T*S0(3), represented as the manifold R^ x SO (3). Using Eqs. ( |5.2|) 
and ( |5.3| ), one can get the following expression for the symplectic two-form uj on the initial 
phase space: 

u = -d{{S X ei)dei)) = -d{eijkSjEikdEii). (5.4) 



It is useful to introduce the variables Ji = eiS = EijSj. For these variables we have the 
following Poisson brackets: 

{Ji, Ejk} = —eijiEik, {Ji, Jj} = —eijk-Jk- 
Note, that we have the equality 

SiSi JiJi- 

The phase space of the top spin model is obtained from the phase space, described above, 
by introducing two first class constraints 

^ = i(52-p2) = 0, (5.5) 
X = - «: = 0, (5.6) 

where 

p > 0, \k\ < p. 

The described Hamiltonian system can be also obtained starting from the Lagrangian 



-\i{ql - 1) - \2{ql - 1) - AgqiQa, 
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where Ai, A2 and A3 are Lagrange multipliers and ^ 7^ is a real constant. The Hamiltonian 
description of the system given by this Lagrangian results in the set of six second class 
constraints 

V'l = - 1 = 0, ^2 = ql-l = Q, ^3 = 91^2 = 0, 

= QlPl =0, ^5 = ^2^2 =0, ^6 = qiP2 + Q2P1 = 0, 



and two first class constraints ( ^.51 ) and (|5.6|) with p"^ = 6^ + , S = Qi x + q2 x and 



63 = Ci X 62, Ci = q^, 62 = 92- The reduction to the surface of second class constraints 
ipi, . . . , ipQ leads, in the end, to the description of the system by dependent variables Eij and 
Si having the Poisson brackets of the form (|5.2| ) and (|5.3|) . 



Consider now the Dirac quantization of the model. Recall that the matrix E can be 
identified with the corresponding rotation matrix (see Appendix C). Let us parameterize 
the matrix E by the Euler angles, E = E{a, (3,'-/), and use the representation where the 
operators, corresponding to these angles are diagonal. In this representation state vectors 
are functions of the Euler angles, and the operators Si and jj are realized as linear differential 
operators, acting on such functions The quantum analogs of the constraints ^p and x 

turn into the equations for the physical states of the system: 

iS^-p^)<^phys=0, (5.7) 
(J3 - K)%hys = 0. (5.8) 

An arbitrary state vector can be decomposed over the set of the Wigner functions, cor- 
responding to either integer or half-integer spins ||3T[| : 



$(«, 13, 7) = (t)jmkD^mM. f3, 7), (5.9) 

where j = 0, 1, . . ., or j = 1/2, 3/2, . . ., and k,m = —j, — j + 1, . . . ,j. The Wigner functions 
Dl^i^ have the properties 

^^Dl^k = jU + ^)D^mk^ S^D^mk = mD^mk^ J^Dlf^ = kD^^^. 

From the representation of the state vector ( p.9[ ) we see that Eqs. ( [5. 71 ) and ( |5.8D have 
nontrivial solutions only when = j{j + 1), and k = k, for some integer or half-integer 
numbers j and k, such that —j<k<j. In other words we get the following quantization 
condition for the parameters of the model: 

p2=j(j + l), K = k, -j<k<j, 0<2jGZ. 

The corresponding physical state vectors have the form 



m=-j 



Thus, we see that the Dirac quantization of the top spin model leads to a pure integer or 
half-integer spin system. 
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Proceed now to the construction of the reduced phase space of the system. As the 
constraints and x have zero Poisson bracket, we can consider them consecutively. Let us 
start with the constraint %p. 



From the expressions for the Poisson brackets (|5.2|) and (|5.3|) it follows that the group of 
gauge transformations, generated by the constraint ?/^, acts in the initial phase space variables 
as follows: 

e,(r) = e, cos(5r) + \{S x e,) sin(^r) + i-S'(5e,)(l - cos(5r)), (5.10) 



S{t) = S, S = \/S\ (5.11) 

Comparing these relations with Eq. (|C.5|) , we see that the transformations under consider- 
ation have the sense of the rotation by the angle St about the direction of the spin vector. 
In matrix notations we can write 

E{t) = ER^{s, St), s = S/S. (5.12) 

Recall that the initial phase space of the system is diffeomorphic to x SO (3). Let us 
consider this space as a trivial fibre bundle over with the fibre SO (3). Gauge transforma- 
tions ( p. 10 ), ( 5.11| ) act in fibres of this bundle. It is clear that the constraint surface, defined 



by the constraint ip, is a trivial fibre subbundle S^ x S0(3). As SO(3)/SO(2) = 5*^, then 
after the reduction over the action of the gauge group we come to the fibre bundle over 5*^ 
with the fibre 5*^. As it follows from general theory of fibre bundles |3^, this fibre bundle 
is again trivial. Thus the reduced phase space, obtained using only the constraint ip, is the 
direct product S*^ x S*^. To find the expression for the symplectic two-form on this reduced 
space we first introduce convenient coordinates in the initial phase space. 

Let TZ{s, n^) be a rotation that transforms the vector n-^ to the vector directed along the 
vector s. This rotation is undefined for 5 = 0, but it is irrelevant for us, because we are 
finally interested only in the constraint surface where S = p ^ 0. We can choose 7l{s, n^) 
in the form 

7^(s, ng) = 7^(n3, (p)Tl{n2, t9), 

where (p and 'i? are the polar angles of the vector S. We shall also use for the 7^(s, 7x3) the 
notation TZ{ip, i)). From the definition of TZ{s, n^) we have 

5 = 57^(s, 713)713. (5.13) 

Using now Eq. ( C.7 ), we can present the matrix R{s, St) as 

R{s, St) = R{s, ns)R{ns, St)R^{s, 7x3). 



Using this representation in Eq. (|5.12|) , we get the following transformation law for the matrix 
F = ER{s,n3): 

F{t) = FR^{ns,ST). (5.14) 
The matrix F defines a new right basis {/j}, with 

fi = FijUj = n^^is, n-i)ei. (5.15) 
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Let A, n and u be the Euler angles for the matrix F, i.e., 

F = R{n3, \)R{n2, /u)i?(n3, z/). 

From Eq. ( |5.14| ) it follows that these angles are transformed under the action of the gauge 
transformations as follows: 

A(r) = A, /i(r) = /i, z/(r) = u — St. 

It is clear that we can use the quantities S, ip, ^ and A, fi, v as independent coordinates 
in the initial phase space of the system. The symplectic two-form a;, written in terms of 
these coordinates has the form 

u = d{S cos dip — Sdu — S cos fid X). (5.16) 

To show that Eq. (|5.16| ) is valid we start with expression ( ^.41) . From Eqs. (|5.13| ), (|C.12| ) and 
i ^JBi ) it follows that 

u = ^d[57^(s,n3)(n3 x fi)d(n{s,n3)f^)]. 
Using matrix notation we can rewrite this relation in the form 

uj = -^d[Str{M3R^{<f,^)dR{<f,^))] + ^d[StT{M3F^{X,fi,u)dF{X,fi,u))]. (5.17) 

Further, Eqs. ( |(J.9| ) and ( |C.11| ) lead to the equalities 

R^{ip,^)dR{ip,^) = {-sm^Mi + cos^M3)dip + M2di&, 
F^(A, /i, i')dF{X, /i, u) = [— sin fi{cosiyMi — sin Z/M2) + cos fiMs\dX 

+ (- sin z/Mi + cos vM2)dfJL + Msrfz/. 

Now using these expressions in Eq. (|5.171 ), and taking into account Eq. (|C.8|) , we come to 
Eq. ( ^361) . 

It is instructive to compare Eq. ( ^.16| ) with Eq. ( [4.6|) . From Eq. ( 5.16| ) we see that the 
restriction of the symplectic two-form to the constraint surface, defined by the constraint if), 
has the form 

UJ = pd{cos'ddip — cos fbdX). 
Using the definition of the variables Jj and quantities Fj^, we see that 

Ji = SFi3. 

From this relation it is easy to conclude that A and fi can be considered as the spherical 
angles, parameterizing Jj, that actually are the polar angles of the vector S with respect to 
the basis {ej}: 

S = S{cos X sin fiei + sin A sin fie2 + cos fie^). 

It is quite clear that the quantities Si and Ji form a set of dependent coordinates in the 
reduced phase space under consideration. Here we have SiSi = JiJi = p^, that confirms 



19 



the conclusion made above that the reduced phase space in this case is S"^ x S'^. Using now 
the experience gained in the preceding Section, we see that the symplectic two-form can be 
written in terms of these coordinates as 



uj — —-z—^i^ijkSidSj A dSk — eijkJidJj A dJk)- 
zp 



Let us turn our attention to the constraint x- K is easy to get convinced that the 
transformations of the gauge group, generated by this constraint act in the initial phase 
space in the following way: 

ej(r) = 7^(e3,r)e^ = ej%(n3,r). 
Sir) = S. 

We see that the gauge group, generated by the constraint acts only in one factor of 
the product 5*^ x 5*^, which is a reduced phase space obtained by us after reduction with 
the help of the constraint ip. Thus we can consider only that factor, which is evidently 
described by the quantities Jj. From such point of view, the constraint surface, defined by 
the constraint is a one dimensional sphere 5*^, where the group of gauge transformations 
acts transitively. Hence, after reduction we get only one point. Thus, the final reduced phase 
space is a two-dimensional sphere S*^. 
In matrix notation we have 

EiT) = R^{n,,T)E, 

that implies 

FiT) = R^{ns,T)F. 
Hence, for the Euler angles A, p, and u we get 

A(r) = A - r, /i(r) = /i, z/(r) = u. 

It is clear from this consideration, that the reduced phase space is described by the indepen- 
dent coordinates ip, {}. The symplectic two-form on the reduced phase space is 

UJ = pd{cosi)dip), 

that can be written in terms of the dependent coordinates Si as 

u = —-^^^ijkSidSj A dSk- 

The reduced phase space we have obtained, coincides with the reduced phase space from 
the preceding Section. Hence the geometric quantization method gives again the quantization 
condition ( [4.16|) for the parameter p, while the parameter k remains unquantized here. 



Therefore, while for this model unlike the previous one, two methods of quantization lead to 
the quantum system, describing either integer or half-integer spin states, nevertheless, the 
corresponding quantum systems are different: the Dirac method gives discrete values for the 
observable J3, whereas the reduced phase space quantization allows it to take any value k, 
such that k"^ < for a system with spin j. 
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Let us note here one interesting property of the system. We can use a combination of 
the Dirac and reduced phase space quantization methods. After the first reduction with the 
help of the constraint ip, the system, described by the spin vector and the 'isospin' vector 
|19| with the components li = —Ji, SiSi = IJi, can be quantized according to Dirac by 



imposing the quantum analog of the constraint x on the state vectors for singling out the 
physical states. In this case we have again the quantization of the parameter k as in the pure 
Dirac quantization method, and, therefore, here the observable J3 can take only integer or 
half-integer value. Hence, in this sense, such a combined method gives the results coinciding 
with the results of the Dirac quantization method. 



6 Discussion and Conclusions 

The plane spin model, considered in Section 3, gives an example of the classical constrained 
system with finite number of the degrees of freedom for which there is no gauge condition, but 
nevertheless, the reduced phase space can be represented as a submanifold of the constraint 
surface. As we have seen, Dirac and reduced phase space quantization methods lead to 
the coinciding physical results for this plane spin model. Moreover, we have revealed an 
interesting analogy in interpretation of the situation with nonexistence of a global gauge 
condition for this simple constrained system with the situation taking place for the non- 
Abelian gauge theories [0 . 

The rotator and top spin models give examples of the classical systems, in which there 
is no global section of the space of gauge orbits. In spite of impossibility to impose gauge 
conditions such systems admit the construction of the reduced phase space. These two 
models demonstrate that the reduced phase space and the Dirac quantization methods (or 
some their combinations), can give essentially different physical results. 

Thus, for Hamiltonian systems with first class constraints we encounter two related prob- 
lems. 

The first problem consists in the choice of a 'correct' quantization method for such sys- 
tems. From the mathematical point of view any quantization leading to a quantum system, 
which has the initial system as its classical limit, should be considered as a correct one, but 
physical reasonings may distinguish different quantization methods. Consider, for example, 
the above mentioned systems. The rotator spin model, quantized according to the Dirac 
method, represents by itself the orbital angular momentum system with additional condi- 
tion singling out the states with a definite eigenvalue of angular momentum operator 
S'^. This eigenvalue, in turn, is defined by the concrete value of the quantized parameter of 
the model: = n{n + 1) > 0. On the other hand, the reduced phase space quantization of 
the model gives either integer or half-integer values for the spin of the system. If we suppose 
that the system under consideration is to describe orbital angular momentum, we must take 
only integer values for the parameter p in the reduced phase space quantization method. 
But in this case we must, nevertheless, conclude, that the reduced phase space quantization 
method of the rotator spin model describes a more general system than the quantum system 
obtained as a result of the Dirac quantization of that classical system. 

The Dirac quantization of the top spin model, or its described combination with the 
reduced phase space quantization gives us a possibility to interpret this system as a system 
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having spin and isospin degrees of freedom (with equal spin and isospin: S = IJi = j(j + l)), 
but in which the isospin degrees of freedom are 'frozen' by means of the condition Iz^phys = 
—k^phys- On the other hand, as we have seen, the reduced space quantization method does 
not allow one to have such interpretation of the system since it allows the variable I3 to 
take any (continuous) value —k restricted only by the condition k"^ < j^, i.e., the operator 
(taking here only one value) cannot be interpreted as a component of the isospin vector 
operator. From this point of view a 'more correct' method of quantization is the Dirac 
quantization method. 

In this respect it is worth to point out that there is a class of physical models, for which 
it is impossible to get the reduced phase space description, and which, therefore, can be 
quantized only by the Dirac method. 

Indeed, there are various pseudoclassical models containing first class nilpotent con- 
straints of the form |2^-[^: 

V' = en-e^„G'^"-^" = o, (6.1) 

where ^j^., are real Grassmann variables with the Poisson brackets 

{^k,^i} = -igki, 

Qki being a real nondegenerate symmetric constant matrix. Here it is supposed that G*^ "*", 
n >2, are some functions of other variables, antisymmetric in their indices, and all the terms 
in a sum have simultaneously either even or odd Grassmann parity. For our considerations it 
is important that constraints { p.l\ ) are the constraints, nonlinear in Grassmann variables, and 
that they have zero projection on the unit of Grassmann algebra. In the simplest example of 
relativistic massless vector particle in (3+l)-dimensional space-time |Q-|22[ the odd part 
of the phase space is described by two Grassmann vectors a = 1, 2, with brackets 

and the corresponding nilpotent first class constraint has the form: 

= ^e,C9"' = 0, (6.2) 

where = diag(— 1, 1, 1, 1). This constraint is the generator of the SO(2)-rotations in the 
'internal isospin' space: 

(^) = cos r + gjj^ (r) = g cos r - sin r. 

The specific property of this transformation is that having ^^(r) and we cannot determine 
the rotation angle r because there is no notion of the inverse element for an odd Grassmann 
variable. Another specific feature of the nilpotent constraint (|6.2| ) is the impossibility to 
introduce any, even local, gauge constraint for it. In fact, we cannot find a gauge constraint 
X such that the Poisson bracket {tp, x} would be invertible. Actually, it is impossible in 
principle to construct the corresponding reduced phase space for such a system. Obviously, 
the same situation arises for the constraint of general form (|6.1|) . It is necessary to note 
here that in the case when the constraint depends on even variables of the total phase 
space (see, e.g., ref. P3); and, therefore, generates also transformations of some of them, we 
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cannot fix the transformation parameter (clioose a point in tlie orbit) from the transformation 
law of those even variables, because the corresponding parameter is present in them with 
a noninvertible factor, nonlinear in Grassmann variables. Therefore, the pseudoclassical 
systems containing the constraints of form ( |6.1|) can be quantized only by the Dirac method, 
that was done in original papers (see also Ref. |]22[ where the BRST quantization 



of such systems was considered). 

Let us come back to the discussion of the revealed difference between two methods of 
quantization, and point out that the second related problem is clearing up the sense of gauge 
degrees of freedom. The difference appearing under the Dirac and reduced phase space 
quantization methods can be understood as the one proceeding from the quantum 'vacuum' 
fluctuations corresponding to the 'frozen' (gauge) degrees of freedom. Though these degrees 
of freedom are 'frozen' by the first class constraints, they reveal themselves through quantum 
fluctuations, and in the Dirac quantization method they cannot be completely 'turned off' 
due to the quantum uncertainty principle. Thus, we can suppose that the gauge degrees of 
freedom serve not simply for 'covariant' description of the system but have 'hidden' physical 
meaning, in some sense similar to the compactified degrees of freedom in the Kaluza-Klein 
theories. If we adopt such a physical point of view, we have to use only the Dirac quantization 
method. Further, the gauge principle cannot be considered then as a pure technical principle. 
From here we arrive also at the conclusion that the Dirac separation of the constraints 
into first and second class constraints is not technical, and nature 'distinguish' these two 
cases as essentially different, since gauge degrees of freedom, corresponding to the first class 
constraints, may reveal themselves at the quantum level (compare with the point of view 
advocated in Ref. p^ ). 



Finally, let us mention also that since the path integral quantization method of con- 
strained systems with first class constraints is based on the introduction of the gauge con- 
ditions for the first class constraints P, ^ , it appears that the usage of this method for the 
systems considered in the present paper is rather problematic. 
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A Pseudoinverse Mappings 

In this appendix and in the next one, we describe the method of dependent coordinates, that 
is a very convenient method to deal with second class constraints (see also ||35||). 



Let A be a linear mapping from an n-dimensional linear space V to an m-dimensional 
linear space W. It is known that for any decompositions 

V = Vo®Vi, W = Wo®Wi, (A.l) 

such that Vq = Kei A, Wi = 1mA, the mapping A induces an isomorphism A' of Vi and Wi. 
Here we have 

A = i^^A'TT^\ (A.2) 
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where tt^^ is the projection of V onto Vi, and is the inclusion mapping of Wi into W. 
Note that decompositions (|A.1|) generate the dual decompositions 

V* = Vq © V*, W* = Wo © W^, 

such that 

WQ=KeTA*, V* = lmA\ 

where A* : W* V* is the adjoint mapping of the mapping A. 
The pseudoinverse mapping of the mapping A is defined as 

= i^^A'-'TT^\ (A.3) 

where tt^^ is the projection of W onto Wi, and is the inclusion mapping of Vi into V. 
From Eqs. (|A.2|) and ( |A.3|) we get 



A-'^A = n, = E 

where 



For the pseudoinverse mapping A^^ we have also 

KeTA-^ = Wo, lmA'^ = Vi, 

furthermore, 



Ker A-^* = Vq*, ImA'^* = W^. 



Choosing different decompositions of the hnear spaces V and W we get different pseudoin- 
verse mappings. Actually, it can be shown that if for a given mapping A we fix decompositions 
( |A.1| ), then there exists a unique mapping A~^, defined by the conditions 

AA-i = S, KerA-^* = V*, 

or by the conditions 

A-^A = Il, KerA-^ = Wo. 

Note that decompositions ( |A.1[ ) can be fixed in the following way. Denote by p the rank of 
the mapping A. Recall that rank A = dimlmA. Let {ea}a=i....,n-p be a basis in Vq = Ker^, 
and {/i"}a=i,...,n-p be a set of elements of V* such that 

In this case we can define Vi as 

Vi = {v e V \ fi^'iv) = 0, a = 1, . . . ,n - p}, 

and we have 

n = 1 - (g) 
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As rank A* = rank A, then dimKer A* = m—p. Let {iy'^}a=i,...,m-p be a basis in Ker A* = 
Wq, and {fa}a=i,...,m~p be a set of elements of W such that 

i^^iU) = 5\. 

The space Wq can be defined as the space generated by the set {fcj}a=i,...,m-p- Indeed, it is 
not difficuh to show that in this case W = Wq © Wi. The mapping S in this case can be 
presented in the form 

S = 1 - ® z/^ 

Treating an arbitrary nxm matrix as the matrix of a hnear mapping from an n-dimen- 
sional hnear space to an m-dimensional hnear space, we associate a hnear mapping with 
any matrix. Hence, we can define the pseudoinverse matrix of a matrix as the matrix of 
the corresponding pseudoinverse operator. We can reformulate the results of the above 
consideration for the case of matrices as follows. 

Let a = \\a'^i\\ he nxm matrix of rank p. Hence, there is a complete set oi n — p 
independent relations of the form 

a^je*Q, = 0, a = 1, . . . ,n — p. 

Since the rank of the transposed matrix equals the rank of the initial one, there is a complete 
set of m — p independent relations 

= 0, r = 1, . . . , m — p. 
It is clear that there are quantities and f^r such that 

IJ, iC 13 = d V rj a = a- 

Define the matrices ||vr*j|| and ||cr''s|| by 

^ j = j - e afJ- j, CF s = s - J s- 

From the discussion above we see that there exists a unique matrix = ||(a~^) V||, defined 
by the conditions 

a\{a~Js = (y\, /i^(a-^)V = 0, 

or by the conditions 

{a-^yy^ = TT^,, (a-i)v/v = 0. 

Consider a special case of skew symmetric nxn matrices. Usually, with such matrices we 
associate linear operators from an n-dimensional V space to the dual space V* . If a matrix 
||ajj|| is skew symmetric then the corresponding mapping A is skew adjoint. A* = —A. If we 



choose in (|A.1| ) 



Wo = V*, W, = V*, 

then the pseudoinverse mapping is skew adjoint and the corresponding pseudoinverse matrix 
is skew symmetric. 
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B Dependent Coordinates 



Let iV be a smooth manifold of the dimension n. Suppose that in some open neighborhood 
U <Z N there is given a set of functions 2;*,i = l,...,m>n such that for any point p & U 
the rank of the set dz^p, i — 1, . . . ,m., is equal to n. In such a case we call the functions 
local dependent coordinates in M. It is clear that if we consider the mapping l : U ^ R"*, 
defined by 

iip) = iz\p),...,z"'ip)), peU, 

then we can consider f/ as a smooth n-dimcnsional submanifold of R™. In fact the case 
when U = M, i.e., the case of globally defined dependent coordinates, is the most interesting 
one. 

Since for any p e U the rank of the set dz\ is equal to n, then there exists a set of 
smooth functions //"j, a — 1, . . . , m — n, such that 

l/^idz^ = 0, a — 1^ . . . ,m — n. 

It is clear that the quantities can be chosen in such a way that the set oi m — n m- 
dimensional vectors ^'^i{p) be linearly independent at any point p & U. 

Let u"-, a = 1, . . . , n, be a set of independent coordinates defined in some open neighbor- 
hood V C U. For the dependent coordinates z'^ we get 

dz' = z'adu", 

where 

dz^ 

z^ — 



The mxn matrix ||-2*a|| is of rank n, and satisfies the relations 

l^"iz\ = 0. 

Let e^a be a set of smooth functions, such that 

1^ ie (3 = p. 

Define the matrix ||7r*j|| by 

Denote by z'^i the matrix elements of the pseudoinverse matrix of the matrix ||2;*a||, defined 
by the relations 

<^ 3 J' 

As it follows from Appendix A, we have the relations 

z^iz'b = z\e\ = 0. (B.l) 

Now we can introduce a set of the vector fields 9/9^*, defined by 

dz' ^ 'du'^' 
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These vector fields satisfy tlie relations 

and they are, in a sense, dual to the differentials dz"^: 

d.-(A)=.V (B.2) 

It is quite clear that the vector fields d/dz^ are independent of the choice of local coordinates 
u*^. An arbitrary vector field X can be represented as 

X = A- A. 

Note that the functions X* in this relation are defined ambiguously. To make them unique 
we subject them to the relations 

/i"iX* = 0. (B.3) 
Analogously, for any one-form = ^pidz'^ the functions (pi are fixed by the relations 

ipie\ = 0. 

Such a fixation allows one to write 

df = ^dz' (B.4) 
dz' ^ ^ 

for any function /, and 

ipiX) = if.X' 

for any one-form Lp and vector field X. 

Suppose now that the manifold is a symplectic manifold with the symplectic two-form 
w, having in the local independent coordinates u"" the representation 

uo = -LUabdu"' A du^. 

Let \\uj'^-'\\ be the matrix composed by the Poisson brackets of the functions z*: 

{z\z^}=uj'^. 

Using the coordinates u°', we can write 

uj'' = uj'''z\zh, (B.5) 

where uj""^ = Hence, we have 

/iW^' = 0. (B.6) 

It is easy to understand that the matrix is of rank n, therefore, the m-dimensional 

vectors form a complete system of the vectors, satisfying Eq. ( |B.6| ). 
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It is clear that the symplectic two-form uj can be represented as 

uj = \jjijdz^ NdzK (B.7) 

The matrix Ht^jjH is supposed to be skew symmetric, and is defined ambiguously. We fix it 
with the help of the conditions 

uji^e\ = 0. (B.8) 

Writing the obvious relation 



UJ. 



a 



and comparing it with Eq. (|B.5|) we conclude that 

Thus, the matrix \\ujij\\ is the pseudoinverse matrix of the matrix Actually, the inverse 

statement is also true. Namely, if we have the representation of the symplectic two-form 
00 of type ( p.7|) , where the skew symmetric matrix \\uJij\\ satisfies relations ( p.8| ), then the 
matrix giving the Poisson brackets of the functions z\ is the pseudoinverse matrix of 

the matrix ||ci;jj||. 

The Hamiltonian vector field Xf, corresponding to the function / on A^, can be written 

as 

.A 



Xf - X\ 



Using Eqs. (jg^), (Q, (|K3|) and we get 



Hence, we have 



and, in particular. 



^ dz^ 



X} = {f,z'}. (B.9) 

Thus, our conventions lead to the same formulae that we have in the case of independent 
coordinates. 

In the present paper we suppose that the one parameter group of canonical transforma- 
tions ^r, generated by a function ip, is defined by the equation 

K{p) = -X4<l>r{p)), peN. 

From this equation for any function / we have 

Kf=f:^{---{f,i^},---,^}- (B.io) 

fc=o ■ ^ ' 

k 
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Consider now the case when the symplectic manifold is a symplectic submanifold of 
an m-dimensional symplectic manifold M. Let Vt be the symplectic two-form on M, and 
L : N ^ M he the inclusion mapping. In this case we have 

UO = L*Q. 

We assume further that the submanifold is a submanifold, defined by the equations 

ip" = 0, a = 1, ... ,171 — n. 
For the functions tjj" we have 

= 0, 

then 

i,*dilj'' = 0. (B.ll) 

Let Z\ i = 1, . . . , m be a set of independent local coordinates in M. It is clear that the 
functions z'^ = i*Z'^ form a set of dependent coordinates in A^. Using the representation 



from Eq. we get 

It is clear that we can put 



dz' = 0. 



dZ^' 

As it was proved in Section 2, at any point p & N the matrix 

is nondegenerate. Using the coordinates Z^, we can write 

n = ^QijdZ' A dZ^, 

hence, 

{Z\Z^} = n'^, 

where \\^'^-'\\ is the inverse matrix of the matrix Thus, we have 

This equality allows us to put 

e\ = (B.12) 

where is the inverse matrix of the matrix 

For the symplectic two-form u we have 

uj = Qijdz^ A dz-' , 
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where flij = L*Qij. Hence, the functions Uij, satisfying Eq. ( [B.8|) , are given by 

For the pseudoinverse matrix we get 

that with the help of Eq. ( [B.12| ) can be written as 

or, equivalent ly, 

{z\ z^} = i*{{z\ p} - {z\ r z^})- 

From this relation for any two functions, F and G, on M we get 

{i*F, i*G} = i* ({F, G} - {F, G}) . 

Thus if we define the Dirac bracket of the functions F and G by the usual relation 

{F, GY = {F, G} - {F, r G}, 

then we have 

{l*F,l*G} = L*{{F,G}*). 

C Three— Dimensional Rotation Group 

Consider the space as a three-dimensional linear space with the canonical basis 






ni = U , 712 = i , ^3 



The scalar product of two vectors u = UiUi and v = ViUi is defined by 

uv = uiVi + U2V2 + Mat's, 

whereas their vector product is 

ux V = eijkUiUjVk, 

where eijk are the components of the totally antisymmetric tensor with £123 = +1. 

The three-dimensional rotation group consists of all nondegenerate linear transformations 
of preserving its orientation and scalar product of the vectors. This group is denoted 
SO (3), and its elements are called rotations. We specify an orientation in R^ considering the 
canonical basis as a right one. Note that the canonical basis is orthonormal. By definition, 
a rotation TZ transforms the basis {rii} to the right orthonormal basis {e-i}: 

Bi = TZrii. 
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We associate the matrix R with the rotation TZ through the relation 



(C.l) 



As the rotation TZ preserves the scalar product, the matrix R satisfies the relation 

RijRkj = Sik, (C-2) 
that can be written in matrix notations as 



RR^ 



1. 



(C.3) 



From this relation it follows that det R = ±1. Remembering that TZ preserves orientation, 
we see that actually 

det R=l. (C.4) 

It is clear that any matrix R, satisfying relations (|C.3D and (|C.4|) , generates a rotation TZ, 
defined by Eq. ([01) . 

A simple geometric consideration shows that the rotation TZ{ni, (f) by the angle (f about 
the direction of Ui is described by the formulas 

7^(r^l,v9)nl = ni, 

TZ{ni, ip)n2 = n2 cos ip + sin ip, 

TZ{ni, ip)n3 = —n2 simp + cosip. 



Hence the matrix R{ni, ip) has the form 

R{ni,(p) = 



1 

COS (/? — sin ip 
sin {p cos ip 



Analogously, for R{ri2,ip) and R{ri^,ip) we get 



R{n2,(p) 



cos ip sin y9 

1 
— sin cos ip 



R{n3,ip) 



( cos ip — sin if 
sin ip cos ip 

V 1 



The rotation by the angle ip about the direction of a unit vector u is described by 

r' = TZ{u, ip)r = rcosip + {u x r) simp + u{ur){l — cosy?). (C.5) 

Hence, we have 

Rij{u, ip) = 6ij cos ip + eikjUk sin ip + UiUj{l — cos ip). (C.6) 

Let Q be the matrix, corresponding to an arbitrary rotation Q. From ( |C.6| ) we can easily 
get 



QR{u,^)Q^ = R{Qu,^). 



(C.7) 



It is quite clear that the Lie algebra of the group SO (3) is formed by the matrices A, 
satisfying the relation 



A + A 



T 



0. 
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A convenient basis in this Lie algebra is formed by the matrices 

dR{ni, (f) 



1/3=0 



dip 

From Eq. ( |U.6| ) we get 
and, therefore, 

[Mi,M,] = ei,kMk. 

A simple calculation gives 

tr {MiMj) = -25ij. (C.8) 
It is clear also that the following equations are valid: 

^^^^^^ = M,i?(n„ ip) = Rim, <^)M,. (C.9) 

dip 

Writing relation (|C4| ) in the form 

^ijk-Ril-RjmRkn ^Imni (C.IO) 

and using Eq. ( |C2| ), we get 

Ril^ijk-Rkn Rjm^lmn- 

This relation in matrix notation has the form 

R^MiR = RijMj. (C.ll) 

Analogously, from Eq. ( |C.10| ) for an arbitrary rotation TZ and arbitrary vectors u, v we get 
the relation 

(nu) X (Uv) = n{u X v). (C.12) 

An arbitrary rotation TZ can be parameterized by the Euler angles with the help of the 
relation 

n{a, (3, 7) = 7^(n3, a)7^(n2, 13)n{n^, 7), 

where < a < 27r, < /9 < vr, and < 7 < 27r. 

Let {ej} be an arbitrary orthonormal basis in R'^. In this case we have 

eiCj = 5ij. (C.13) 

Expanding the vectors ej over the canonical basis {rij}: 

Bi = E.jTij, (C.14) 

we get the matrix E = {Eij). From Eq. ( |C13D it follows that 

EE'^ = 1. (C.15) 
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The basis {cj} is a right one if and only if 



det E = 1. 



(C.16) 



Hence, any right orthonormal basis in is connected with an element of the group SO (3). 
On the other hand, any matrix E, satisfying relations (|C.15|) and (|C.16|) , specifies a right or- 
thonormal basis in by Eq. (|C.14| ). Thus, we see that there is a one-to-one correspondence 
between the set of all right orthonormal bases in R'^ and the group SO (3). 
Note that condition (|C16D can be written as 



~ijki 



or as 
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